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\^ ', In many physical situations the behavior of a quantum system is affected by 

^ ■ interaction with a larger environment. We develop, using the method of influence 

QhI functional, how to deduce the density matrix of the quantum system incorporating 

Qh! the effect of environment. After introducing characterization of the environment by 

1^ • spectral weight, we flrst devise schemes to approximate the spectral weight, and then 

a perturbation method in field theory models, in order to approximately describe the 
environment. All of these approximate models may be classified as extended Ohmic 
models of dissipation whose differences are in the high frequency part. 

The quantum system we deal with in the present work is a general class of har- 
monic oscillators with arbitrary time dependent frequency. The late time behavior of 
the system is well described by an approximation that employs a localized friction in 
the dissipative part of the correlation function appearing in the influence functional. 
The density matrix of the quantum system is then determined in terms of a single 
classical solution obtained with the time dependent frequency. With this one can 
compute the entropy, the energy distribution function, and other physical quantities 
of the system in a closed form. 

Specific application is made to the case of periodically varying frequency. This 
dynamical system has a remarkable property when the environmental interaction is 
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switched off: Effect of tfie parametric resonance gives rise to an exponential growtfi 
of the populated number in higher excitation levels, or particle production in field 
theory models. The effect of the environment is investigated for this dynamical 
system and it is demonstrated that there exists a critical strength of the friction 

for the parametric effect. In the model of periodically oscillating field coupled to 
a system quantum field, it is verified that the parametric effect occurs in medium, 
with a somewhat diminished rate, if the relaxation time scale of the system field 
towards thermalization given by the friction term is larger than the time scale of 
the coherent parametric amplification. The effect persists until the back reaction 
against the periodic oscillation stops particle production. The resulting energy dis- 
tribution of produced particles described by a universal function deviates from the 
thermal one, having an average energy that exponentially increases with time. The 
dynamical system driven by the parametric oscillator thus maintains and does not 
lose its quantum nature even in thermal bath. In the present work analytic formulas 
of how physical quantities behave at late times both in the high and in the low tem- 
perature regions are given, along with results of numerical computation displaying 
time evolution. 



2 



1 Introduction 



Any system of physical interest has interaction with the rest of a much larger 
system, and one must clarify effect of the environment on the system under study. 
The problem can be bypassed in some idealized situations in which the effect of 
the environment is considered to be too small and essentially negligible. In the 
other unique case one deals with the entire universe and in this case there is no 
clear separation between the system and the environment. Except in these special 
cases one often has a clear idea of how to separate the small system and the large 
environment, with some kind of interaction between them. The important issue one 
deals with is then to estimate the effect of environment without knowing too much 
of its property, usually difficult to precisely specify. Even more crucially, in many 
physical problems effects of the environment are of central importance. As examples 
of such, and only those we come to recognize immediately, we may list Brownian 
motion, quantum tunneling and decoherence phenomena. 

The first purpose of this paper is to develop a general formalism in a special and 
important subclass of the system-environment interaction: a quantum system of 
time dependent harmonic oscillators q interacting with the environment variable Qa 
via bilinear coupling: q J2a CaQa ■ There must exist infinitely many environmental 
variables distinguished here by the index a for the environment to behave as a large 
reservoir. There already exists a large body of literature that focuses on this class 
of models, but we introduce a spectral characterization of the environment system. 
In this way we can devise useful schemes of approximation for the environment 
consistently with general principles. 

Extension to an infinite set of time dependent harmonic oscillators q^ for the 
system in question has a variety of applications in many diverse topics. This class 
of models, for instance, includes an interesting case of field theory models in which 
the system field ip interacts with the environment field x via the Yukawa type of 
interaction Lagrangian density: >Cint = ^ f , with fj, a dimensional coupling 
constant. The bihnear composite field plays the role of the environmental variable 
Qa in this model. The index k for the system in question is taken here to be the wave 

— * 

vector k of mode if we are interested in a homogeneous environment. We propose 
to study this type of the system-environment interaction beyond the often-discussed 
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simple model of the quadratic type, oc (p Ca Xa ■ We have some interesting results 
for this new class of field theory models at one loop level. 

The second purpose of this paper, and actually the original one that motivated 
us towards this problem, is to examine effects of thermal bath on particle production 
under periodic perturbation. There are a few cases that particle production under 
periodically oscillating field may become important issues in cosmology. It was 
recognized in cosmology that periodic coherent field oscillation may give a major 
source of entropy production in our universe, perhaps the most obvious candidate 
being inflaton oscillation right after inflation. This specific cosmological problem is 
called the reheating problem The essential part of the reheating problem belongs 
to the class of models of time dependent harmonic oscillators set up for the first 
purpose. The time dependent frequency in this case is given by a sinusoidal function 
oscillating like cos(mg t) . The inverse period is the inflaton mass in the reheating 
problem. What has been discovered in recent study p| — [0 of the reheating and 
related problems is a very rapid and efficient mechanism of particle production via the 
parametric resonance. The effect drastically changes the old picture |T^ of particle 
production based on the simple inflaton decay. However, in the usual discussion 
of the parametric effect one ignores the presence of the cosmic medium. Here are 
some important, unsolved problems: without the environmental interaction there is 
no randomness or entropy among produced particles. Another application of the 
parametric resonance in cosmology is decay of very weakly interacting field such 
as the Polonyi field, which has a danger of creating too much entropy when they 
decay at late times. In these problems it is not immediately clear how the effect of 
the environment affects the rate of particle production, or even the presence of the 
parametric effect itself. We shall considerably advance our understanding on these 
problems in the present work. In particular, we prove for the first time that the 
parametric effect can occur in medium, and give a critical condition for this. 

We must stress in the beginning limitation of our approach: that our formalism is 
based on the possibility of the clear cut separation of the system and the environment. 
Hence it is difficult for us to discuss the important phase of the reheating problem: 
namely when the environment itself rapidly changes by created particles, ultimately 
dominating the rest of the environment. We thus cannot discuss the entire history of 
reheating. Nevertheless, we may give a definite picture of how particles are produced, 
including their energy distribution at the instant of their creation prior to interaction 
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among created particles themselves. 

On the other hand, our formalism is flexible enough to cope with many states of 
the environment. For simplicity and as our initial step of investigation, we treat the 
environment in thermal equilibrium characterized by a finite temperature T = 1/(3. 
Needless to say, this includes the quantum ground state as a limit of T ^ 0. 

The rest of this paper is organized as follows. In section II we first review briefly 



the influence functional method of Feynman and Vernon [^. We discuss how the 
correlation function a(r — s) that appears in the influence functional is identified 
to the real-time thermal Green's function in the operator formalism. We then im- 
plement this method by introducing the Lehmann spectral weight for the case of 
bilinear coupling between the system and the environment variables. The spectral 
weight contains information of both the coupling and excitation spectrum of the 
environment, hence we call this the response weight. The correlation function in 
the influence functional is completely characterized by this response weight. We 
organized this section from the point of the inluence functional, to stress that this 
method, when supplemented by the spectral representation, becomes a powerful tool 
to characterize the system-environment interaction. 

In section III we devise, consistent with some general properties of the response 
weight, approximate models of the weight consisting of simple poles of relevant an- 
alytic functions. This model extends the Ohmic dissipation model often used in the 
literature and has the response weight of the form, r{uj) = ^uj /(^) . tj is the friction 
constant and the cutoff function f{x) ^ as a; ^ oo. A general argument is given 
to support the Ohmic behavior at low frequencies: r{uj) (x uj as ^ 0. The second 
class of approximations we discuss is perturbation in terms of the coupling /i in the 
field theory model given by Ci-^t above. A full response weight at one loop level is 
given in this field theory model. If the main dissipation of the system field (p occurs 
via X + ^ X, as so happens typically at high temperatures, the field theory model 
may also be regarded as an extended Ohmic model. The cutoff function differs much 
in various models discussed in this work. We compute in all of these approximations 
the friction constant that becomes important in subsequent sections. 

In section IV we introduce, as a model of the quantum system, time dependent 
harmonic oscillators. If one solves this case for generic time dependence in a closed 
form, it should be possible to deal with much more general class of models for the 
system. We explain in detail a crucial assumption for the correlation: the localized 
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friction approximation for the dissipative part of the correlation function Q;/(r) ~ 
r]6'{T). This approximation is excellent for discussion of the late time behavior. 
Under this approximation the reduced density matrix of the quantum system is 
completely determined in terms of a single classical solution u{t) of the system, 
obeying il +uj'^[t)u = , with the initial condition, u{0) = , u{0) = 1 . t^^(t) here 
is modified from the original one in a known way. We explicitly compute the density 
matrix, the energy distribution function, and the entropy of the system, taking as 
the environment the one in thermal equilibrium. The energy distribution function, 
when written in terms of a scaled energy divided by the average energy, is given 
by a universal function that has a single spectral shape parameter. We also check 
how all this works in the special case of time independent harmonic oscillators, and 
compare our result with previous works |]T^ — [0. Some comment on the early 
time behavior is also made. 

In section V we discuss effect of thermal bath on the quantum system of frequency 
CO subjected to periodic oscillation, C,{t) = C,o cos(mgt) . Without the environmental 
interaction this system exhibits the exponential growth of population in oscillator 
excitation levels when the parameters of the oscillation are in one of the infinitely 
many bands of the parameter space {uo/m^ ^^o/'^^)- We clarify the condition under 
which the parametric resonance enhancement persists despite interaction with the 
thermal bath: If the time scale in which a disturbance of the system variable is 
relaxed to equilibrium is larger than the growth time scale of level population, then 
the exponential growth occurs with a somewhat diminished rate compared to that 
in vacuum. The condition of the exponential growth also gives a critical strength of 
the friction above which the parametric effect does not occur. For quantum fields 
coupled to a periodic field oscillation this condition limits the momentum range of 
particles produced by the parametric resonance. Furthermore, although the energy 
distribution of produced particles deviates from the thermal one, it has a long tail 
in the high energy side. The average energy is shown to exponentially increase with 
time. We present analytic formulas in the high and in the low temperature regions, as 
well as results of detailed numerical computation including time evolution of physical 
quantities. 

Section VI gives a summary of the present work. Some technical detail related 
to numerical computation is relegated to Appendix. 
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2 General framework: influence functional and 
spectral representation 



The Lagrangian of our problem consists of three parts: using the system variable 
q and the environment variable Q, it is 

L = L,[q]+U^,[q,Q]+LQ[Q]. (2.1) 

We take for the system-environment interaction the bilinear term: 

Lint = -qY^ CaQa ■ (2.2) 
a 

In later discussions we mainly focus on the time dependent harmonic oscillator for the 
system dynamics, but for the time being, namely in the present and the next sections, 
we may choose any Lq for the system dynamics. Since we do not specify the nature 
of the environment variable except that Qa^ are independent, the bilinear coupling 
Lint above may be regarded as a general one. For instance, the environment variable 
Qa may be a composite operator made of more fundamental fields, as discussed in the 
field theory model in the next section. In any case it is important to have infinitely 
many indices a coupled to the single system variable g, because the environment 
must behave as a large system compared to the g— system. 

It is very useful to introduce the influence functional, following Feynman and 
Vernon |rT|. The basic idea is that one is interested in the behavior of the g— system 



alone and traces out the environment variable altogether. We define here the influ- 
ence functional by convoluting with the initial state of the environment. To do so 
we assume that initially we may take an environment state uncorrelated with the 
system. In the path integral framework the influence functional is thus obtained 
after integrating out the environment variables: 

^[g(r),g'(r)] ^ j VQ{r) J VQ'{t) J dQ, J dQ[ J dQf J dQ) 5{Q j - Q)) 
K ( g(r) , Q(r) ) K* ( q\r) , Q\r) ) p,(Q, , Q[) , (2.3) 
K{q{T),Q{r)) = exp ( «^o[g] + «^int[g , Q] ) , (2.4) 

5o[Q] + Sint[g , g] = f dT{ Lq[Q] + Lint[g , Q] ) • (2.5) 

JO 

It is a functional of the entire path of the system g(r) and its conjugate ^'(t). 
Pi{Qi , Q'i) is the initial density matrix of the environment, which can be any mixture 
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of pure quantum states. What deserves to be stressed is that one does not observe 
the final state of the environment, hence integration with respect to the final values 
oi Qf = Q'jr is performed here. 

Once the influence functional is known, one may compute the transition ampli- 
tude and any physical quantities of the g— system by convoluting dynamics of the 
system under study. For instance, the transition amplitude is given by 

I Vq{T) I Vq'ir) J dq, J dq[ J dqj J dq'f 

r{qf) rii) :F[q{T) , g'(r) ] e^^^-^l^'J ^(g,) ^(g^) , (2.6) 

where ■^'s are wave functions of the initial and the final q'— states and S[q\ is the 
action of the gr— system. 

The form of the infiuence functional is dictated by general principles such as prob- 
ability conservation and causality. Feynman and Vernon found a closed quadratic 
form consistent with these. 



•^[?(T),g'(r) 



exp 



f dr r ds ( ^{T)aR{T - s)^{s) + i ^{T)ai{T - s)X{s) ) 
Jo Jo 



(2.7) 



with e(T) = q(r) - q'(r) , X{t) = q{r) + q'{r) . (2.8) 



Thus two real functions Q;i(r) are all we need to characterize the system-environment 
interaction. These are defined here in the range of r > . The fact that ctj depends 
on the difference of time variables, r — s, is due to the assumed stationarity of 
the environment. Without this stationarity these may depend on r , s individually: 
ai{T,s). 

It would be useful to recall that if one takes as the environment an infinite 
collection of independent harmonic oscillators in thermal equilibrium of temperature 
T= then 

Mt) = Ec>£H^)> «/(^) = Eca4"^W> (2-9) 



a 



(«)/' ^ _ ^-v,/'^'^«^ cos (a;ar) 



i?^(r) = coth(^)^-^, (2.10) 
(a) sin(a;„T) ... 

with sign(T) the signature function, sign(T) = 1 (r > 0) , — —1 {r < 0) . We 
extended the r range to — oo < r < oo such that the function in r < is consistent 
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with eq.( p.l5| ) later discussed; a\"'\—T) = a;-"''(r) . If one Fourier-transforms 

a(«)(r) = 4")(r) + ta'f\T) , (2.12) 

with 

/oo 
dr a^^^T) e''^\ (2.13) 
-oo 

then one finds the famihar form, 

Emergence of the thermal propagator is by no means an accident. By compar- 
ing the path integral formula and the operator formula, one can indeed derive the 
equivalence of to the real-time thermal Green's function; 

a('^)(r)=tr (p^T[g„(r) g,(0) ] ) . (2.15) 



Pj3 is the thermal density matrix e~'^'^'5/(tr ^^"^ ) of the environment. In order 
to avoid confusion, we explicitly denoted the environment operators by the tilde 
notation. The key point for this identity is the form of the bilinear coupling, Ljnt = 
— q J2a CaQa ■ Taking as the initial environment an arbitrary quantum state, not 
necessarily limited to the thermal state, gives the infiuence functional of the form, 

J'[q{T),q'{T)] = J VQ{t) J VQ'{t) J dQ, J dQ,' p,{Q, ,Q/) 

J dQj J dQf'6{Qf~Qf')exp[t{So[Q]-So[Q']) 

-I f drJ^cA e(r) {Q + Q')a{r) + X(r) {Q - Q'Ur) 
2 Jo V 



(2.16) 



which is to be equated to 

exp[- [ dr [ ds^{T) {an^T ,s)^{s) + iai{T ,s)X{s))] . 
Jo Jo 



Here 



SolQ] = I drLgiQ) (2.17) 
Jo 

is the action of the environment. Functional derivatives with respect to ^(r) and 
X(r) yield various combination of expectation values of the operators, Q{t) ± Q'{t) . 
A suitable combination of these leads to the required identity ( ^.15| ). 
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It should be evident that the fundamental relation (|2.15|) holds for any initial 
Q— state if one takes the relevant density matrix pj. Equivalence of the path integral 
formula to the corresponding operator formula is used to extract useful characteri- 
zation of the environment, as will be explained. 

We now introduce the Lehmann spectral representation. Without explicitly solv- 
ing interaction with the environment, this spectral knowledge is adequate to predict 
the behavior of the system under study. Respecting some general properties of the 
spectral function given below, one may also devise useful approximate schemes for 
the spectral density. This approach is thus very useful to characterize the environ- 
ment without knowing too much of it. 

Assuming the thermal density matrix for the initial Q— system, one obtains, by 
using QaiT) = e*^'3'^ Qa(0) e~*^'3'^ and inserting a complete set of eigenstates \n) 
of the Hamiltonian Hq, the thermal real-time Green's function: 



(2.18) 



zG',(r) =tr iT[Qair)QaiO)]pf3 



El e 



n ,m 



The Fourier transform defined by 



/TT 



yields with ujmn = - E„ 



dTGJT)e'' 



r) .(2.19) 



(2.20) 



IT E l(^IQa(0)|m) 



1 

2^ 



+ ie - bJmn 
1 



duj' 



iU -ie + LUmn 

ie 



1 



1 Taiuj') 



^ 1 - e-/^'^ uj -io' + ie 1 - e'^'^ 
(l_e-^-) Y: \{n\QM\m)\'e-^^-6{u;-u;^^^) 



UJ — UJ' — It 

1 



n,m 



El e 



-fSEi 



The quantity 

a 

we call the response weight, in this paper. 



(2.21) 
(2.22) 
(2.23) 

(2.24) 
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We shall discuss general properties of the response weight r{uj) shortly. It would 
be instructive to recall beforehand for a single harmonic oscillator for Q, 

ra{u) = ^ {5{u~uJa)-S{u + Ua)) . (2.25) 

ZUJa 

This can be computed by using Qa{0) = (a^ + al) in terms of the creation 
and the annihilation operators of harmonic oscillator and summing over the energy 
levels, En = {n + ^) Ua ■ Despite the initial appearance of the temperature factor 
P it disappears in the final form of the response weight, as expected from previous 
results. For a collection of harmonic oscillators in thermal equilibrium, the response 
weight is thus 

r{u;) K^' - ^l) sign (^) = E ( ^(^ " ^a) - S{u; + (2.26) 

The response weight in this form was introduced in [l^, but we stress that here it 



is defined in the most general way, using the spectral representation for arbitrary 
environment. 

The response weight r{uj) thus defined characterizes the response of elementary 
excitation in thermal environment that may couple to the system variable. Thus 
only this quantity needs to be specified to determine the behavior of the g— system. 
By definition it obeys a positivity constraint, 

r(tu)sign(cj) > 0. (2.27) 

Furthermore we often assume as an acceptable system of environment that the high 
frequency limit of the response weight is bounded such that 

r{Lo) 0[u''] ascj oo , (2.28) 

with e > 0. If this condition is not satisfied, one must use instead of r{uj) a better 
behaved r{uj)/uj'^ with appropriate n > 1 in subsequent dispersion integrals, and 
must introduce as a subtraction term a polynomial in u of order < n — 1. It is 
important to note that the response weight is an imaginary part of some analytic 
functions. With the usual type of retarded and advanced Green's functions, 

'^('^) = rd:^G«M + 13557 G-^l-)' (2-29) 

G''A^^) =r '•'"'> =r'^ , , (2.30) 

j-oo 27r uj — Lu' ±ie Jo vr co^ — uj' ± ietu 
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one has 

r{uj) = 2QG^{uj) = -2'^G^{uj) . (2.31) 

The retarded (advanced) ) is analytic in the upper (lower) half cu— plane. 

A sum rule on the spectral function is usually taken granted. This holds if Qa is 
a canonical variable, and takes the form of the normalization integral, 



duj uj ra{uj) . (2.32) 

We do not consider this constraint any further, because our response weight ( p.24{ ) 
contains the unknown coupling and an unspecified number of dynamical variables. 

Although the real-time formalism is a natural consequence of the influence func- 
tional method, it is sometimes more convenient to use the imaginary-time formalism. 
If r{uj) obeys the convergence condition, eq.( |2.28| ), then there exists a unique analytic 
function G{z) of a complex variable 

(2.33) 

J-oo ZTl Z — Uj' 

that gives the real-time thermal Green's function G^'^iu) for real u as limiting 
functions: 

G^(cu) = G{uj + ie) , G^{uj) = G{uj - ie) . (2.34) 
The imaginary-time thermal Green's function G defined at discrete Un is simply 

GM = GM , (2.35) 



with LUn = i^TT n/ (3 , n = ±1 , ±2 , ■ ■ ■ . This is a well known result |T^. We 
note that with the uniqueness of the analytic function, the analytically continued 
imaginary-time Green's function gives a convenient means of computing the response 
weight as a discontinuity: r{uj) = (^G{uj — ie) — G{uo + ie)^ /i . We shall use this 
fact to compute the response weight in field theory models. 

With this machinery the correlation function that appears in the Feynman- 
Vernon formula is given by 



2 



oo 



ai{T) = I rfcur(cu)e-*^^sign(r) 

-oo 



/ (iu;r(ct;) sin(u;r) sign (r) , (2.36) 
Jo 



1 I3uj 
(^r{t) = o du coth — r{uj) e~ 

Z J — oo Z 



/ (ico'r(co') coth — cos (cur) . (2.37) 
Jo 2 
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Their Fourier transforms give the famihar dispersion relation at finite temperatures, 

aniuj) = TT coth ^ r{uj) , (2.38) 
2 ^^,,00' tanh ^ aR(u}') 

TT Jo iO^ — iO' ^ 

In the next section we shall discuss two useful approximate schemes for the re- 
sponse weight r(uj\ one based on simple meromorphic forms of this function and the 
other based on perturbation in field theory models. 



3 Models for environment 



It is usually difficult or even impossible to know details of the environment, yet 
one must be able to extract essential features of the reservoir in order to predict the 
behavior of the system under study. Based on the response weight introduced in 
connection with the infiuence functional in the preceding section, we discuss two ap- 
proximate schemes to evaluate this quantity. These two methods are complimentary. 
In one method we follow the spirit of the S-matrix theory: Instead of computing the 
response weight from fundamental microphysics we approximate relevant analytic 
functions in terms of simple pole terms whose parameters may be used to character- 
ize the environment system. In the other method we solve the fundamental problem 
by perturbation theory. What one can compute in this case is necessarily limited, 
for instance only the lowest order computation is readily tractable. Nevertheless one 
may gain some insight of how to relate macro-quantity to microphysics parameters. 
There is also a freedom of how to model the environment. We introduce non-linear 
models of the system-environment interaction in field theories. This type of models 
seem to have a variety of areas of application, some of which we hope to cover in 
subsequent works. 

Ill A. Meromorphic approximation for the response weight 

The critical observation to devise useful scheme of approximation is that the 
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response weight r{uj) defined for uj real is a limit of some analytic function, 

rH= limjsfiz), (3.1) 

where f{z) = — 2 G^{z) is analytic in the upper half complex z— plane. The simplest 
approximation for G^{z) is to replace its complicated analytical structure by a set 
of simple poles consistent with general principles. Moreover, to extract the late 
time behavior of the correlation function a(r) it is sufficient to determine the low 
frequency limit of the response weight r{uj), as will be explained. The minimum set 
of poles thus might be a good approximation. Even the more complicated case of 
singularities such as the branch cut may be described by a continuous distribution 
of poles. 

It turns out that the oddness and the positivity on the positive real axis of uj 
constrains much of pole parameters of the response weight r{uj). Poles thus appear 
in a pair ai z = Wi and W2 where W2 = — wl , with residues q related by C2 = — . 
Using the minimum set of poles, one has 



1 



UJ — W UJ + W UJ — w* u + w* 
4cf27 UJ 



(3.2) 

4 J' + ^^n^ 

with SfJw 7^ 0, and for < , 3fJc > and = . There are three real parameters, 
9fJc , "^w = r2 , = — 7/2 . We call this approximation a quartet model. 
The exception occurs for = 0; 

.m^ ^ ( ^ 1 ^ 



r^^^(a;) 



uj — i'~)/2 uj + i'~)/2 ) 
2auj 

(3.3) 



+ 774 ■ 

We call this case a doublet model. In the literature it is called the Drudes model 



and has been repeatedly discussed |jTj]. We shall not treat this model in any detail. 

From this response weight the correlation function that appears in the Feynman- 
Vernon formula becomes 

ag) {to) = TT coth ^ r» {uj) , (3.4) 

(4)rr) = -27rce-^l"l/=' sin(fi|r|), af (r) = - vr a e'^l'l/^, (3.5) 



a 



00 



"iJV) = 4cfi7 i coth ^ — — ,^^2^^ , ^ cos(^r) , (3.6) 



2 (cj2 _ ^2 + + Q2^: 

anHr) = 2a [ duj coth — — ■ cos(co'r) . (3.7) 

JO 2 + 7 /4 
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To gain more physical insight, it is useful to explore the late time behavior at 
l/Vt (with f2 7) for the quartet model. As r — > 00, af\T) behaves as rj5'{T) 

with 

r] = -2 r dTTaf\T). (3.8) 
Jo 

The leading asymptotic behavior of the correlation function ^/(r) in the quartet 
model is thus given by 

„«M~,y(r), „ = (3.9) 

Putting these forms of a;/(r) into the Feynman- Vernon formula (p.7|) , one gets the 
friction term of the form, 

Vqir)-^q{T), (3.10) 

ignoring a surface term. We shall treat the surface term rigorously in later applica- 
tions. 

The quartet model actually contains the frequency shift term — jj^q^ij) as a sub- 
leading asymptotic term, which is difficult to extract unambiguously. The frequency 
shift turns out to correspond to the mass renormalization in field theory, and is not 
an observable quantity. On the other hand, the localized friction term above is very 
important, and with the localized approximation for Q;/(r) above the real part of 
correlation is given by 



a 



00 



g^(r) = / duj coth(— )r(^)(cj) cos(cjr) , (3.11) 
JO 2 

r(^)(o;) ~a;r'(^)(0) = ^o;. (3.12) 

TT 

This approximation for the response weight r^uj) is often adequate in determining the 
behavior at asymptotic late times. But the frequency integral on uj is not convergent 
unless one uses the full form of r(u;) such as the one given in the quartet model. We 
shall discuss the high frequency cutoff later. 

The form of the response weight, r{uj) oc cu, is often called the Ohmic model of 
dissipation in the literature. In this and the next section we are giving our micro- 
scopic foundation of this form of the behavior modified at some high frequency cutoff 
Vt. Although the Ohmic model does not come out of fundamental principles, its low 
frequency behavior appears in a variety of interesting models. To show this we first 
consider a class of the response weight that comes out of the harmonic oscillator 
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bath; 

rH = ^DH. (3.13) 



We replaced the discrete sum in the previous formula eq. (|2.26| ) by a continuous 



integral introducing the density of states: D{uj). In many systems of physical inter- 
est both in condensed matter and particle physics the density of states behaves as 
D{uj) oc cu^ at C(j ~ 0. With c{uj) — > a constant as — 0, this means the Ohmic 
behavior, r{uj) oc u. In analytic approach to the response weight one may consider 
approximate models in which poles are continuously distributed; 

rico) = ^co r dnfi^). (3.14) 

We assume the convergence of the Q integral. The cutoff function /(^) can be taken 
either the one in the quartet model or in the doublet model. The Ohmic behavior 
r(ci;) oc at cij ~ is then a consequence of the presence of a gap; fio > . 

Departure from the Ohmic behavior appears at some frequency scale and if this 
scale also coincides with the cutoff scale of the high frequency contribution, one may 
regard the model as an extended Ohmic model. On the other hand, models that 
contain several distinct frequency scales may considerably differ from the Ohmic 
model. Our main focus in the response weight in the present paper is on the Ohmic 
model modified at high frequencies: extended Ohmic models. 

Ill B. Field theory model and perturbative analysis 

The infinitely many environment variable Qa can actually be a composite field 
of more fundamental fields. In order to discuss a fully relativistic field theory, we 
also consider the case in which the system variable qk is a Fourier component of 
some scalar field (p{x) put in a homogeneous medium. As a simplest field theory 
model of non-linear system-environment interaction, we thus take a Yukawa type of 
interaction with a new bose field x{^)'^ 

Lint = d^X Cint , C,nt = " ^ ' (3-15) 

with fi a parameter of mass dimension. V is the normalization volume. By identifying 
the Fourier-mode with 

Qk e''-' = -^=i= ( a, + aU ) e^^^ , (3.16) 
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we deduce the environment variable coupled to it, 

E CaQa = ^ / d'xx\x) e-''-\ (3.17) 

Henceforth we regard the field x &s a fundamental variable of the environment, 
assumed to be in thermal equilibrium. 

The correlation function a(r), or more conveniently its Fourier transform a{uj), 
can be computed as in the real time formalism at finite temperatures. To lowest 
order in fj,"^ 



/ d'xe'^-Hi [t[x\x,t)x\^.s)]pp) (3.18) 



2 

is calculable in terms of the x~propagator in the momentum space; 

1 2m 



uo'^-P- ml + it e/^-^fc - 1 

X. 



K^'-uji)]. (3.19) 



This is of course identical to the Fourier transform of the one we are already familiar 
with, aR(r) + zQ;/(r). We omit the suffix x for the x~niass in what follows. 

As in the T = field theory it is easiest to compute the discontinuity instead of 
the full diagramatic contribution. Physically the most transparent way to compute 
the imaginary part of the self-energy diagram, 53 , is to use the analytically 
continued expression [l^ from the imaginary-time formalism, as outlined in the 



preceding section. The analytically continued imaginary part, or more precisely the 
discontinuity, directly gives the response weight; 

r{ijj) = 2^Ii{uj) = 2u}T{u}) , ^Ii{u}) = ^, {Y[{uj - it) -Ii{uo + it)) , (3.20) 

where T{ijj) is interpreted as a decay rate in thermal medium. 

In the subthreshold region of < one loop contribution is given by 

'^Yi{uj) = r dE {n{E) -n{E + uj)) , (3.21) 



ijj k / 4m^ 

where 

^(E) = (3.23) 
is the Planck distribution function. With 

n{E) - n{E + to) = n{E){ 1 + n{E + to)) - n{E + uj){l + n{E) ) , (3.24) 
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the imaginary part ( p.21| ) for |ci;| < A; is a sum of the two contributions, X + X 



and its inverse processes that occur in thermal medium. The factor 1 + n represents 
the effect of stimulated boson emission. Near a; = 0, 

yU^ a; 



^H~3^T^(^*)' ^* = - VF + 4m2 . (3.25) 



This in turn gives the friction coefficient, 



V = nr'{0) = -^niE,). (3.26) 



As emphasized in |18], the decay rate 



r(0) = ^r'(0) = |- (3.27) 

can be interpreted as the relaxation rate of thermal disturbance of the y?— field when 
departure from equilibrium is introduced. At high and low temperatures the friction 
thus computed behaves as 

"-i^TTO-'^--'' .-ge-v-^/.-.,(r^O). (3.28) 



On the other hand, for on > vP + 4m^ relevant physical processes are ip ^ 
X + X- Since 

( 1 + n{E) ){l + n{E + uj)) - n{E)n{E + cu) = 1 + n{E) + n{E + cu) , (3.29) 
the imaginary part of the self-energy diagram is given by 

(a;) = P ( 1 + 2n{E) ) , (3.30) 



as UJ ^ oo. When divided by oj, the first term here r(ci;) = 9'n(cc;)/a; is ~ 
/i^/(327ru;) , which is the decay rate for — > XX vacuum, including the effect of 
prolonged lifetime at c<j ^ m^p . The response weight given by the imaginary part of 
the self-energy at one loop is shown in fig.l for a set of the parameters ( , m , T ) , 
both in the sub- and super-threshold regions. Also shown here is an approximation 



to this weight, eq.(|3.34]) given below. 



Ill C. High frequency cutoff and extended Ohmic model 
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The truncated form of the response weight is often adequate in discussion of most 
physical quantities, especially at high temperatures. But the frequency integral for 
the correlation ani^r) is not convergent if r{uj) oc uj all the way up to a; = oo. 
Actually there exists a high frequency cutoff; 

nVl^ 1 

r^^\uj) -, for quartet model , (3.32) 

,,2 

r^^Ucu) for field theory model, (3.33) 

Svr 

assuming ri ^ 7 for the quartet model. We subtracted the T = contribution in 
writing the limiting form r^'^^u) of the field theory model. 

The cutoff frequency appears quite different in the two cases, at Q in the quartet 
model and at T in the field theory model at finite temperatures. Moreover, the high 
frequency tail is also different and we observe a very fast convergence of the temper- 
ature dependent term in the field theory model. In the literature the straightforward 
cutoff at some a; = is also used frequently |jl3[, [|1^]. It is important, however, to 
note that there may be temperature dependence of the cutoff Q. We believe that 
in some realistic applications the relation Q oc T holds at high temperatures. For 
instance, if one computes excitation energy in the form of mass correction in finite 
temperature field theory, one usually finds that Q ^ gT at high temperatures, with 
g some dimensionless coupling constant. In the rest of discussion we assume this 
proportionality at high temperatures, whenever we need realistic estimates. At very 
low temperatures this relation does not hold, and the cutoff Q reflects some intrinsic 
scale of the environment, to which we have to pay attention. 

It appears that the constant weight at a large u, r{uj) ~ /i^/(167r) , in the field 
theory model must be retained at low temperatures. This constant term does not 
contribute to the imaginary part of the correlation ^/(r) (see eq.( |2.3(j| )). We may thus 
use the subtracted weight, r{uj) — r(oo) for ^/(t), which behaves asymptotically 
like 0[l/c<j^] as cij — * 00. But if the high frequency contribution of the response 
weight with r{uj) /i^/(167r) dominates over the low frequency contribution in the 
uu integral of the real part of the correlation, it does not fit into the approximation 
scheme of the localized friction, at least in the present simplest form. If we do not 
deal with this case and remain in the realm of the localized friction approximation, 
the field theory model belongs to an extended Ohmic model with a cutoff Q ~ 



k"^ + 4m^ . Physically this means that we identify the main process of dissipation 
in the field theory model as X ~^ f ^ X instead of (f ^ x ~^ X ■ 



19 



A better approximation of the field theory model is to replace the one loop 
response weight by a continuous function that well reproduces the correct behavior 
both in the a; = and in the a; = oo regions. If one ignores a mismatch in the 
threshold region, u — k ^ A;^ + 4m^ , then an example is given by 

= ^ \^ eoth ( -^^T^ \ , (3.34) 

^ ' 167r 2w + A^o A; \^ 4T j ' ^ ' 

No = e^^'+^/(2T)^;L_ (335^ 

As illustrated in fig.l, this gives a reasonable description of the response weight 
except in the threshold region. The threshold region alone may be described by a 
quartet type of model, but we shall not elaborate on this problem any further. 

Thus in all models discussed in this section the response weight has the Ohmic 
behavior at low frequencies. Differences in the models are in the form of the cutoff 
function /(^) and the cutoff scale f2; f{x) oc in the quartet model, f{x) oc l/x"^ 
in the doublet model, and f{x) oc 1/x in the one-loop field theory model. If one 
focuses on the late time behavior, the form of the cutoff function is not critical unless 
the environment temperature T <^ Q . In subsequent applications we use a simpler 
cutoff of the form, e~^^^ , whenever necessary in analytic estimates. On the other 
hand, for numerical computation we cutoff various uj integrals at ^2 as a simple means 
of computation. This should be sufficient practically when application is limited to 
high temperatures. Comparison with other forms of the cutoff has been attempted 
in a few sample cases. We first checked that results of numerical computation shown 
below are insensitive to the value of the straightforward cutoff scale Q beyond some 
value, Q > 50 in our unit. Use of the quartet model that has the same friction 
constant and the same effective maximum value of the response weight as in the case 
of the cutoff method neither did produce any recognizable change of our results. 



4 Time dependent harmonic oscillator in thermal 
bath 

We have so far discussed how to characterize the environment using the infiuence 
functional. We now apply this formalism to the dynamical system governed by the 
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Lagrangian of time dependent harmonic oscillator; 

L,^\e-\u.\t)ci\ (4.1) 

The main purpose to treat this class of models is to clarify effect of the environment 
on the parametric particle production. Besides this direct goal this class of models 
can cope with many interesting cases of application, some of which we hope to cover 
in subsequent works. We therefore explain how to derive basic formulas somewhat 
in detail. 

The basic quantity we focus on is the path integral over the system variable g; 
Vq{T) j W(r)exp(z%(r)]-^%'(r)])^[g(r),g'(T)] 
= TV / Vi{r) j VX{t) exp [\ [ dr (^(t) X{t) - uj\t) i{T)X{T) ) 

- f dr r ds ( e(r)ai^(T - s)i{s) + ze(r)a/(r - s)X{s) ) 1 , (4.2) 
^0 Jo 

with ^(r) = q{'T) — q'{T) , X{t) — q{T) + q'{T) . The effect of the correlation a{T — s) 
introduces non-local interactions among the system variable q, in general difficult to 
deal with. Even diagrammatic expansion may not work with the presence of time 
dependence in u!^{t). But the most important cases can be discussed using the 
approximate localized friction for ^/(t — s), to which we shall turn. In this case one 
can explicitly perform the path integral in terms of a classical solution for ^(t). The 
problem is then reduced to analysis of a rather simple classical system. 

IV A. Localized friction 

The notion of localized friction arises by focusing on the late time behavior of 
q;/(t). By the late time here we mean the asymptotic time that is much larger 
than any time scale of excitation in the environment; t ^ 1/ (maximal excitation 
frequency) 1/Q. The limiting form is then given by 

/■oo 

aj{T) = — / (ia;r(cc;) sin(tc;r) r]5'{T), (4.3) 

JO 

/•oo 

drTai{r)^nr'{0), for r(a; = oo) = , (4.4) 

Jo 

OiR{r) - / da;c^/(— ) coth(— ) cos(a;r) . (4.5) 

The essence of this approximation is the low frequency truncation to the response 
weight, 

r(u;) =u;r'(0) = -a;, (4.6) 

TT 
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such that f{x) fa 1 . We already computed in the preceding section the friction 
constant 77 in a variety of models, specifically in the quartet model and in the field 
theory model. 

Although the dissipation given by a/(r) is local by the time scale of l/Q, the 
damping time scale of the noise kernel aR^r) may differ. It can be shown that at 
high temperatures of T > f2 the locality of the correlation holds for r > giving 

«il(T)~^5(T). (4.7) 

TT 

On the other hand, at low temperatures of T < ^2 the local approximation may be 
inadequate and a high frequency cutoff is needed. As already discussed, we replace 
in this case the response weight by the more appropriate one; for instance in the 
quartet model 

r(a;) = ^a;/(^), /(^) = (^^^TTWiWtW^ ' ^'"^^ 
or by the simple frequency cutoff at uj = Q with f{x) = e^^, as is often practiced. 

With a/(r) oc S'{t), the exponent factor in the path integral becomes, including 
the mass renormalization effect. 



2 



Xi]\ -'-n^fXfsign{t) - 1^ dr £ dsar)aR{r - 3)^8) ■ (4.9) 



+ 2 

A remarkable feature of this formula is that a part of the action excluding the 
term is a local integral described in terms of a renormalized frequency wr{t) and 
the friction term oc rj. This much is enough to considerably simplify calculation of 
the transition amplitude of the g— system. 

Path integral over the system variable q is standard; it can be decomposed into 
contributions from the classical path and from a Gaussian fluctuation around this 
classical path. Let us first write down the equation for the classical path of ^{r), 
which is obtained by functional differentiation with respect to X{t); 

This equation is best analyzed by introducing a new ^(t), 

y(r) ^ e^^/' ^r) , (4.11) 
(^2+-Rir)-^-^)yir) = 0. (4.12) 
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It is convenient to use two types of solution u{t) and v{t) with different boundary 
conditions; 

w(0) = 0, v{t)^0. (4.13) 
Actually if u{t) is known, the other is given by 

with c a constant to be determined by normalization. 

The normalization of these solutions is fixed by specifying the value of the con- 
stant Wronskian, 

u{r)v{r) - v{t)u{t) = - u{Q)v{Q) = u{t)v{t) . (4.15) 
For instance, by setting this to the initial shifted frequency 



a;o = y a;o' - J , = ^fl(O) , (4.16) 

one may have, consistently with the Wronskian relation, 

m(0)=c^o, v(0) = 1. (4.17) 

In the rest of this paper we leave formulas as general as possible without specifying 
the normalization condition, but in places where we specifically need normalization 
constants we shall use this normalization. 
The classical ^ path is then 

«='M = «^^-""-'"^^«-$|^""^- 

Fortunately, we do not need details of the other solution Xc\{t); only the boundary 
values Xi ,Xf are needed for the X part. 

The classical action is thus computed after getting , ^/ from ^ci{t), 

^s., = -^ef-^^^-w^,^f-'-vCfXf 

- -xJ^d^ + "^)+^f^e-^^A , (4.19) 
2 ' \^'^2 v{0) ' u{t) y ' ^ ^ 
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V 



u 



2 [ dr [ 

Jo Jo 

2 f dr I 
Jo Jo 

/* dr r i 
Jo Jo 




r) (2t-T-s)/2 



(4.20) 



(4.21) 



W 



ds 



u{t)v{s) + u{s)v{t) 



ttR(r-s)e-''(*-"-^)/^ (4.22) 



The second step of computation on the Gaussian fluctuation around the classical 
path is trivial for SX{t) = X{t) — Xci(r) , because the action is linear in this 
variable. This simply gives a delta-functional containing a linear term in 5^(t); 



Since the boundary condition is S^{0) = 6^(t) = , the only solution that makes 
the argument vanish is the trivial one, S^{t) = 0. This makes the rest of SC,{t) 
integration also trivial. The net result is then a determinant factor, 



The Feynman- Vernon formula (|2.7|) respects unitarity to arbitrary orders of the 
correlation Q;(r). This makes it possible to handle the determinant factor above by 
a total normalization factor N{t), which is fixed in the end of calculation by the 
condition of the density matrix, tr p^^^ = 1 . 

Finally, we note that although non-local dissipation is beyond the scope of the 
present work, there are physical situations that inevitably go beyond the local friction 
approximation. In the extended Ohmic model that has a physical high frequency 
cutoff fl, the short time behavior at r < l/fi cannot be trusted, because the picture 
of the local friction given by a;/(r — s) oc 6'{t — s) fails. If the time scale of the 
system change is ujq, it must be that ujq < Q. Otherwise what takes place in the time 
scale of T < 1/ujq cannot be trusted. The non-local dissipation becomes important 
in the case of uq > Q . 

IV B. Reduced density matrix, energy distribution, and entropy 

A useful quantity for determining the behavior of the quantum system under the 
environmental interaction is the reduced density matrix p^^^ obtained by convoluting 
with the initial state of the g— system; 




(4.23) 




(4.24) 



P 



'^^\Xf,^f)=N{t) / V^r) / VX{t) / d^, / dX, 
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■po{ii.X,) eyLY>{iS[q{T)]-iSy{T)])T[q{T),q\T)\ 



(4.25) 



We take for the density matrix of the initial state, 



. .coth(,.„.„/2) exp[-^coth(gf^)g-^ta.h(g|^)Xn. (4.26) 

which gives a thermal state at temperature To = V/^o- Wc use this initial ^— state as 
a convenient means of computation, since we can treat both the pure ground state 
(in the limit of Tq 0) and a thermal state. 

ghtforward Gaussian integration then gives the reduced density matrix. 



A straigl 

A 



B = 



C = 



1 /«(0).2 r.-l 

^ + ^coth(^)(^)^e-*- 

2 AujQ ^ 2 



it)- 



2K 



K = S + e 



L = C + 



1 m(0) L I ( 7] u{t) 

2^K^ 2 \2~ Mtj ) ' 

^*^coth(^)fl + i,(!^ + 44)^ 



^ ^2 + ' 



'A,B,C) = / dcu coth(^)r(o;)(a,6,c)(a;), 
^0 2 



b{u) = 
c{oj) = 3fJ 



70 W 



/o w 



(0) 



M(i) Jo 



The normalization is determined from probability conservation, 

/oo 
cigp(«)(2g,0) = l, 
-oo 

leading to 



(4.27) 
(4.28) 
(4.29) 
(4.30) 
(4.31) 
(4.32) 
(4.33) 
(4.34) 

(4.35) 
^(*-)^ .(4.36) 

(4.37) 



lA 



p^^^ {Xj,if)=2\\-eM-^Xj-Bij + tCXfij]. 



(4.38) 
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This is a fundamental result of the present work: The reduced density matrix is 
determined by two factors: a single classical solution u{t) (and related v{t)) for the 
system, and the response weight r{u) for the system-environment interaction. The 
only approximation we need to derive this formula is the local form of a/(r — s). 
The diagonal part of the reduced density matrix with q = q', namely with ^ = 0, 

is 

p(-f^) (Xf = 2g , 0) = 2 J- exp[ -AAq^]. (4.39) 

V TT 

The width of the Gaussian peak of 1/ VA is a measure of how the system behaves, 
but we know more details in the form of the complete density matrix. 

It is sometimes useful to transform the density matrix in the configuration space 
to the Wigner function fw{x ,p), 

fw(x,p) = l^^d^pi^){q^x+^-,q'^x-^)e-^P^, (4.40) 

Expectation value of the number operator, namely the occupation number, in terms 
of the reference frequency, taken to be the initial one ujq, is calculated most easily 
from the Wigner function; 

which consists, except the trivial |, of two terms, the term ljq/{16A) from the 
Gaussian width of the diagonal density matrix element and the rest from the kinetic 
term — We shall discuss the rate of particle production in more detail when we 
turn to specific applications, namely the two cases of constant and periodic frequency 

One can work out the frequency distribution of populated harmonic oscillator 
levels. Suppose that one asks the probability of finding the system in the unit 
energy interval around u — ^ + x'^ . We get the energy distribution function f(uj) 
from the Wigner function; 

^^^^ " /^^^^^'^)'^^^-y-T^') (^-^^^ 

= — V»exp -— — + — ^ , 4.44 
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^ = (— + 



4^ ac+uji ^ 4^y^' 



Here 



1 

TT Jo 



■ u> cos 



is the modified Bessel function, having 

1 



as w ^ oo) , 



1 (as — > 0) . 



(4.45) 



(4.46) 



(4.47) 



V 27ri/; 

This distribution differs from the thermai one and has some interesting features. 
At a; — >• oo, it decreases as 



A 



D = 



AA AC+ujI ( AA AC + ujI^^ AA\ 



1/2 



+ 



B J 



(4.48) 
, (4.49) 



assuming 2; 7^ 0. The high frequency tail of this distribution decreases more rapidly 
than that of the Boltzmann distribution, e~^'^. On the other hand, at a; = 



The average energy of the distribution is computed as 

roo 

(u) = / dujujfiuj) 
Jo 

UJqB ( AA 4^2+^.2 



(4.50) 



AA V <^o 



AuoB 



This average satisfies the relation. 



(4.51) 



(4.52) 



as expected. 

The exceptional case of 2; = with an exact exponential form is actually im- 
portant, as will be realized in thermal equilibrium. This is possible only when 
C — 0, and V AB — ojq/A , giving 



with T ^2B, 



(4.53) 
(4.54) 
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The normalized spectrum shape written in terms of a scaled energy, x — u)/ (a;) , 
is characterized by a single parameter 6; 



Mx;5) = -^e-^/^'loi^^^^x) (4.55) 



e"^/-^' r de expf-x^--^ — cos^l , (4.56) 
Jo 0^ 



with which the energy distribution is given by 



S ^ .J^ hA^^Sll^y . (4.58) 

This function has the following limiting behavior: 

Mx;6) e"^ as5 ^ 1, (4.59) 
fo{x;5) ^S(^ ^ 0. (4.60) 

The first moment (x) happens to be the same in the two limits of 5 = 1 and 5 = 0. 
The main difference of these two limiting distributions is in the low x part. In fig. 2 
this universal function is plotted for a few choices of S. We call S the spectral shape 
parameter. 

An important measure to explore the behavior of the quantum state under the 
environment action is the effective entropy of the system under study. We define as 
usual the entropy of the system as (henceforth we omit the index R for the reduced 
density matrix) 

,5 = -trplnp, (4.61) 

where the trace operation is performed on the system variable q. A useful device to 
compute the logarithmic matrix is to convert it to a power series according to the 
formula like 

The arbitrary power tr is calculable, 

/n 
n dQi p(^) {qi + q2,qi-q2)--- P^^^ {Qu + qi,qn-qi) 
i 

{Vb + Va)^-{Vb-Va)^' ^' ' 
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A more efficient way to compute the entropy is ffist to observe, with the presence 
of the power series expansion, independence of the factor C in the entropy formula. 
This justifies neglect of the C term (by setting C = 0) in the density matrix ( [4.38| ) 



for computation of the entropy, and only for this purpose. The next step is to 
identify this density matrix with C = as an equivalent harmonic oscillator system 
of frequency u) under a thermal bath of temperature 



The equivalence is possible only for B > A and is established by the parameter 
relation, 

u = aVAB, tanh^ = W^. (4.65) 
With this identification the equivalent temperature and the entropy are given by 



T = — 7 — 1= 1= ^ , (4.66) 

S = -J^-\n{l-e-(^n 

= -(x-l)ln--^ + ln^, (4.67) 



with 



x=^^ (4.68) 

Clearly the effective frequency u) differs from the average {uo) computed from the 
distribution /(cu), eq. ( [4.511 ). 

When A '^B holds, considerable simplification follows; 



AVAB, (4.69) 
f ~ 2i3, (4.70) 
S ~ iln^ + l. (4.71) 

Comparison with the energy distribution function f{uj) makes it clear that the 
thermal equilibrium state is possible only in the limit of 

C = 0, ^-0, (4.72) 
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in which T 2B. The second condition imphes that 

^ = 2^/1^0: (4.73) 

The energy quantum huo must be much smaller than the thermal energy ksT. 

When B = A, the density matrix is factorizable into a product of tp{q) ip*{q') , 
since tr p*^ = 1. This is a pure quantum state, having the entropy S* = . 

Existence of the equivalent harmonic oscillator basis for i3 > ^ actually implies 
much more than a convenient means of the entropy computation. First, ii B < 
A, tr p" > 1 for all n > 2, which means that despite an apparent probability 
conservation tr p = 1, some eigenvalues of p are negative; violation of the positivity 
of the density matrix. Secondly, we have found the basis of the states that diagonalize 
the density matrix at any instant of time: It is the Fock basis \n) oi the equivalent 
harmonic oscillator of frequency u) with diagonal elements; 

It is thus important to check in any physical applications the inequality, 

i3 > ^ > , (4.75) 

which is necessary for the positivity of the density matrix p. In all practical applica- 
tions we shall discuss subsequently, the required inequality is obeyed, but when the 
picture of the local dissipative kernel is doomed to fail, this inequality, B > A\s also 
violated. These occur either in the very short time behavior at r < , or in the 
case of the low frequency cutoff, Vt < ujq. This has been checked in numerical compu- 
tation of the harmonic oscillator model mentioned in the next subsection. Thus the 
model of the local friction fails very badly violating the positivity condition, when it 
does fail. One has to be very cautious about the validity of the local approximation. 

Finally, let us mention for comparison results that comes out when interaction 
with the environment is switched off: riuj) = or = . The form of the density 
matrix (|4.38|) is the same as before with 

B = coth^(^) A=-^ coth(^) , (4.76) 
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When the system is in a pure quantum state initially with /5o = oo, it always remains 
in a pure state subsequently with A = B . Thus the entropy S vanishes, while particle 
production can occur; (A^) > . 



IV C. Time independent harmonic oscillator and early time behavior 

A special case belonging to our class of dynamical systems is the independent 
collection of simple harmonic oscillators in which ujj^^s are time independent. Al- 
though the simplest, this type of models is a fundamental model of the system- 
environment interaction and used in various physical situations: to name just two 
of them, quantum Brownian motion [0 — |T^. and decoherence problem [|T^. We 



shall recapitulate main computational results in our way of approach, to the extent 
that comparison with application in the next section is useful. In the end of this 
subsection we discuss, using these results for time independent uij^, the early time 
behavior for the general case of time dependent C(j^(r). 

The basic classical solutions in the case of a constant frequency are, with uj^ = 

'^o = ^0 - T ' 

sin ( cun (t — r) ) , „ x 

uir) = sin r) , v{t) = \ °1 . . 4.78 

sm^cuo t) 

Most of the relevant quantities can analytically be computed, except the u integral 
involving the factor coth ^ at finite T. The late time and/or the high temperature 
limits are fairly simple, and to avoid too much complexity of formulas we only record 
the high temperature and the late time behavior. In these limits we may use 

S^^g, (4.80) 

,^c^m_i(v_m]. (4.81) 

2B u{t) 2 V 2 u{t) I ^ ' 

What we mean here by high temperature actually covers a wide range of inter- 
esting region. What should be compared with the temperature T here is the energy 
scale of the system, ujq. The other quantity rj is taken to be 0[g'^T] with g some 
coupling < 1. Thus as far as T ^ cuq is satisfied, one can regard the system in the 
high temperature region. 
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A straightforward computation gives in the T ^ oo hmit; 



lJ2 rp 

A ^ B ^ -, C ^ 0. (4.82) 

ST ' 2 ' ^ ' 

With these, physical quantities behave as 

{N) ^ —, 5^1n— + 1, (4.83) 

(u) ^ T, f{u) ^ (4.84) 

The entropy obeys the law, e"^ ~ (A^) . Furthermore the density matrix given by ^ , i3 
and C here agrees with the thermal density matrix of the environment temperature 
T in the T — oo limit. 

These results indicate that the final state of the g— system is insensitive to its 
initial state, as seen by independence of Pq, and it is the universal state at the same 
temperature as the environment. The approach of (uj) and 5* to this asymptotia is 
fast, with the exponential e~^^ factor in correction terms, as illustrated in fig. 3 and 
fig.4 for a few choices of the friction rj and the initial system temperature Tq. In these 
and all other subsequent numerical computations we set the cutoff scale = 50. All 
other frequency and energy scales are determined relative to this unit. The final 
physical quantities are also insensitive to how one treats the high frequency cutoff 
and even to the values of t] and Tq. The asymptotic values in these figures do agree 
with the analytic formula, eq.( [4.83| ). These are presumably an expected result, but 
computations here show the power and the simplicity of our method. 

We finally discuss the behavior of physical quantities at early times, going back 
to the general case of time dependent frequency uj]^(t). Strictly speaking, our for- 
malism only allows a precise estimate of the behavior of the g— system at late times. 
Nevertheless it would be of great interest to extrapolate results to earlier times. This 
is indeed the way people discuss the early time behavior in the literature As 
noted previously, the very early time behavior at r < 1/Q cannot be trusted in 
the approximation of the localized friction. This means that to us the early time 
behavior refers to 

^ < r « - . (4.85) 

To this discussion we do not need to solve the classical equation for u{t) for the 
case of the time dependent frequency. Only this function for constant cj^, eq.( [4.78D 
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is needed, since time dependence appears in higher powers of t for small t. Thus the 
leading terms as i — > 0+ are 



A, B, C^O[f], 
u{0) 1 v{0) 



1 



(4.86) 



(4.87) 



u{t) t ' v{0) 



t 



Using these, we find that 




77 
2 



(4.88) 



Physical quantities behave as 




(4.89) 




-ita„h(«^)expl-§^(l + |l)ta„h(/^)l 



(4.90) 




(4.91) 



Without the system-environment interaction {1] = 0) the quantum system exhibits 
both vacuum and thermal fluctuation, while the entropy of the system is given by 
the thermal contribution alone. 

Thus the system starts at its given temperature Tq = I/Pq. Correction to this 
leading term arises only in powers of time t, except the rj dependent terms in all 
quantities excluding the entropy. These r] dependent terms appear to show a energy 
flow from the environment without accompanying entropy. It persists even at Tq = 0, 
the case of a pure quantum state. 

In the case of time independent harmonic oscillators the system finally approaches 
the same thermal state of the environment smoothly. We show in fig.3 and fig.4 result 
of numerical computations for constant cuj^, to indicate how the occupation number 
(fig.3) and the entropy (fig.4) evolve with time. 
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5 Quantum system under periodic perturbation 



We finally specialize to the case in which a quantum system field ip is coupled to 
a coherent field oscillation S,{t) oc cos(mgt) . (For a lack of notation we use the same, 
already used Greek letter ^ to denote the external oscillator field in this section.) 
The coherent field may be viewed as an aggregate of zero-momentum particles with 
some kind of precise coherence, is thus the mass of these bosons. In our present 
work this field oscillation is regarded as given, hence we do not discuss field damping 
due to particle production of the <y9— field. Our main concern here is in the effect of 
thermal bath on the particle production caused by ^—oscillation. 

Under this circumstance the harmonic oscillator variable is qk = ^-^^^y ( + 
a}_^, ) , where the creation and the annihilation operators are those of Fourier k— 
mode (oc e*'^'^) of the system field '^[x). For the oscillator-system interaction we 
take the quartic coupling given by 

^/eV, (5.1) 

primarily because this case has been analyzed in considerable detail [0 without 
taking into account the environment effect. The /c— mode variable then obeys the 
evolution equation, 

{ ^ 



+ u^r) gfc(r) = , uj\t) = P + + g^^^ cos^m^r) . (5.2) 



^0 is the amphtude of oscillation. The standard form of this type of equation is called 
the Mathieu equation and is usually written in dimensionless units; 

\ dz 



^+h + 2e cos(22) q{z) = , (5.3) 



z = mi:T, h = ^ + 2^, e = ^. (5.4) 

4m| 

It is well known that this quantum system exhibits instability in infinitely many 
band regions of the parameters of ( \k\ ,gC,o) , or {h,6) . It has been realized 



that this instability gives rise to particle production ||^. The most important recent 
development is that for large amplitude oscillation particle production and associated 
field decay is greatly expedited, typically ending in a time scale of O[10 — 100] times 
1/m^ However, in all initial investigations so far no systematic estimate 
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of the environmental effect has been attempted. This is precisely what we wish to 
do in the present section. 

Taking over the formalism developed in the previous sections, we ask how the 
classical solution u{t) and v{t) evolve with time r, from which we can compute all 
physical quantities of the quantum system. Let us first recall that these functions 
satisfy the Mathieu equation modified by the frequency shift and the friction: 

[-^2+^"(^) + ^^"-j) ^ir) = 0. (5.5) 

The frequency shift is in proportion to the cutoff: 600"^ oc Q , hence is divergent in 
field theory models. As usual in renormalizable field theories it is more appropriate 
to regard the renormalized mass + 6uj'^ as a physical parameter. This function 
u{t) is related to the original mode function g(r) by 

n(r) = e''"/^ g(^) . (5.6) 

Limiting our analysis to modes within the instability band of the modified Mathieu 
equation, we note that two linearly independent solutions are either growing or 
decaying according to 

e^"^e-p(r), e" ^'"^^ i?(r) , (5.7) 

where A > is the dimensionless growth rate and both P{t) and R{t) are periodic 
with period of 27r/m^. Needed solutions must satisfy the boundary condition, m(0) = 
, v{t) = , hence these are linear combinations of the growing and the decaying 
solutions. With the normalization given in eq. (|4.17|) , leading asymptotic behaviors 
of these are 

u(t) ~ e^""^' Pit) , vit) = - !^(^^)!^ ^ e~ ^'"^^ R(t) , (5.8) 

u{t) 

with both P{t) and R{t) bounded functions. 

Note that the growth rate A here is related to the solution of the modified Mathieu 
equation, eq. (|5.5| ), and not of the original one ( ^.2|) . Thus we should keep in mind 
that the growth rate A does depend on the friction rj: \ = \{ri) . Suppose for instance 
that the renormalized y9— mass vanishes: + 600"^ = 0. The t] term in eq.( ^.5|) 
then has an effect of lowering the band level, since the h parameter in the standard 
Mathieu equation is given by 

h = h-^ = ^ + 2e-^ (5.9) 
4m| m| 4m| 
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in this case. Although it is in general difficult to estimate the influence of the friction 
f] on the growth rate X{ri), the general trend is obvious: The lowering effect of the 
band level tends to increase the growth rate. In the extreme case of a very large rj 
the h parameter may become negative, and the parameter is clearly in the instability 
region. Needless to say, in perturbation theory we later work with, the leading A 
term is independent of 1]. Of course, the real criterion of the parametric amplification 
is derived only after we work out time evolution of physical quantities. 

The exponential growth of the original mode function g(r) is possible for 
2X{ri) > 1]. We observe two competitive factors for the growth: the rate of the 
parametric amplification 2X{ri) against the friction 77. As will be shown in subse- 
quent computations of physical quantities, the friction, when it is small, does act to 
diminish the parametric particle production, but does not wipe out the parametric 
effect. Thus, instead of the blocking factor, it is more appropriate to view the friction 
as the inverse time scale for the system to be driven towards thermalization. Hence 
if the friction is small enough, the parametric amplification never loses against the 
friction. 

Although implicit, it is important to remember v{t) depends on t through the 
boundary condition. Thus it is more appropriate to write this as v{t ; t). An example 
of these classical solutions, u{t) and v{t) normalized at r = t and r = respectively, 
are shown in fig. 5. 

In subsequent discussion we present both analytic and numerical computation of 
time evolution of physical quantities. In analytical computation we can give only 
the limiting behavior. On the other hand, detailed numerical computation has been 
performed, which covers a variety of situations. Some of these results are shown in 
figures. In analytic estimate we concentrate on the asymptotic late time behavior 
for which r ^ max. {l/rj , l/(Am^) ) . With the asymptotic behavior of the classical 
solution, ■u(r) and v{t), the three quantities that appear in the reduced density 
matrix ( 4.38|) of the (^—system are 

2K u{t) 2 1, 2 u(t) ) ■ ' ' 

The quantities K ,A ,C are already defined in eq. ( [4.31| — ^4.36| ) and have asymptotic 
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behaviors, 

A = 0[1], C < 0[Max. (e-^™«* ,e"''*/2)] , = 0[Max. ( e-^^^"?* , e"^*) ] . 

(5.13) 

Thus, we conclude that A ,B ,C behave asymptotically as 

^ = 0[Min. (l,e-2^'"«*+'^*)], B = 0[1], C = 0[1] . (5.14) 

We omitted dimensional scales of these quantities such as Uq and r], since these differ 
from case to case, depending on the parameter range of T , ujq , , r] . For the 
reference frequency we take Uq = ujji{0) : the initial value of the frequency. We also 
note that these asymptotic forms are sensitive to the initial state, dependence of 
To = l/Po, only via the K term of eq. ( |4.31| ). 

The formulas so far are valid for any response weight r(u;). We now specialize 
to the environment system that may be approximated by the localized friction with 
r{uj) = ^uj X (high tu— cutoff). In particular, we are very much interested in 
the high and low temperature limits of various quantities. As in the preceding 
case of constant frequency the high temperature here effectively means that T ^ 
Max. {uo ,1] ,Q) . 

First, in the high temperature limit of T — * oo the frequency integral is well 
convergent and one has 

/■°° , 1 / / NX 27r ^, , , ^ ^, 

/ auj coth -— uj cos{uj{t — s)) — o{T — s) , (5.15) 
Jo 2 p 

Thus in terms of powers of T, A = 0[1/T] , B = 0[T] , C = 0[1] . Let us first 
focus on the case of strong parametric resonance for which 2Am^ ^ t] such that the 
effect of parametric resonance is well operative. In this case 

r]T Jo ^ u{t) ' 



dTe~"^'~^U \ [ ) ! ) / dr e-'>^'-^U^Y \ , (5.17) 
^ ^u{t)v{Q)' j \Jo \{0)' J 
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Note that the expression for C is different from the constant uj^ case discussed at the 
end of the previous section, because in the present case u{t) exponentially grows. The 
integrals here are crudely estimated by noting that u{t) increases rapidly at the time 
scale of l/(Am^), much shorter than l/rj, starting from r = 0, while v{t) decreases 
at the same time scale. We checked this behavior numerically by integrating the 
Mathieu equation. Thus 

From these asymptotic forms we find at high temperatures 

(.) ,5.22) 

1 p-W(2H) 

f\u) ^ ^ — (5.24) 



Both the number {N) and the average energy {u) increase exponentially with time, 
the rate being somewhat diminished from the value taken without environment. Note 
also that a rather non-trivial enhancement factor oc T appears in the pref actor. 

We show in fig. 6 — fig. 8 time evolution of various quantities resulting from nu- 
merical computation at relatively high temperatures in the cases of both the small 
and the large friction. The adopted parameters relevant to the original Mathieu 
equation are h = 1 , 29 = 0.75 , and correspondingly k ~ 0.5 , g^Q ~ 1.2 . 
This belongs to the intermediate ^—amplitude region giving A ~ 0.185 in vac- 
uum. The global growth rate computed from the modified Mathieu equation is 
2A(?7) — ~ 0.17 in the case oi r] = 0.2, which fits well with numerical results of 
time evolving physical quantities. Observe in these figures that the case of a large 
friction given here hy 7] = 1.0 does not yield the parametric amplification, as inde- 
pendently checked by 2A(?7) — rjr^— 0.69 < . An outline of how we numerically 
computed A,B ,C is described in Appendix. An example of time evolution of these 
basic quantities is shown in fig. 9. 

The energy distribution ( ^.24|) exaggerates the low energy part, because it does 
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not yield the correct cu ^ limit. A better approximation is 

1 p-'^/{2M) 

fiu) = (5.25) 



7r4C2 + ^2 Am^ V 2 u{t) 

The low energy cutoff coc <^ (cc;) (the average). Although the energy distribution 
differs from the thermal one, the dominant component is in the high energy side. 

In this model the total system combined with the environment may gain energy 
only from the oscillating field ^{t). The energy input rate from the coherent oscillator 
is estimated by 

It = " IM sin(2m^i) . (5.27) 

The rate itself contains oscillating factors, but as shown in fig. 10, the integrated 
energy ^ ^ 

{E,,,) = - dr sin(2m^T) , (5.28) 

always increases after transient times of t < 0[7r/(2m^)]. A deep reason how 
this comes about is the cooperative interference of the two oscillating components, 
sin(2m^i) and the one in 1/A(x{^)^. If u{t) monotonically increases without 
oscillation unlike in this case of the parametric resonance, there is no net gain of the 
energy {E^ot) from the oscillator. With regard to this it is interesting to note that in 
the case of a large friction the quantity 1/A does not exponentially grow, but has a 
regulated oscillatory behavior, which results in a linear, instead of the exponential, 
growth of the energy flow (£'tot) as illustrated in the inlet of fig. 10. 

In a sense the quantum system of the parametric amplifier maintains in one way or 
another its quantum nature even at high temperatures, provided that 2X{'r)) m^ — r] > 
. This is in sharp contrast to the case of simple harmonic oscillator put in thermal 
oscillator bath, for which we have seen that the late time behavior is the same 
thermal state of the environment. It might be possible to experimentally observe a 
quantum behavior of the parametric amplifier at high temperatures. 

Despite the suppression of the growth rate from 2Am^ to 2Am^ — r], the environ- 
ment interaction does not erase the parametric effect if 

2X{r])m^>r]. (5.29) 
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The derivation of this growth condition constitutes one of the most important con- 
clusions of the present work. Let us look into the meaning of this condition that 
guarantees the exponential rate of particle production in thermal bath. Since rj is the 
relaxation rate of (p— field disturbance towards thermalization, the condition simply 
implies that the thermal relaxation never catches up the parametric amplification for 
a sufficiently small friction. At the same time the condition implies that there exists 
a critical strength of the friction rj above which the parametric effect does not occur. 
The equation 2X{ri) = rj thus gives the critical strength of the friction constant 
T]c. We wish to stress that although the suppression of the growth rate might have 
be anticipated from a naive consideration , it was not clear in the past how 
one can verify the parametric amplification itself, starting from quantum mechanical 
principles. In our view the significance of the friction t] is in its role as the relaxation 
rather than destruction of the coherence. Hence if the relaxation time scale is larger 
than the amplification time scale, the parametric effect wins over the decoherence 
due to the environment. 

The growth condition eq.( |5.29| ) can be quantified. We shall base our estimate 
of this inequality on the field theory model that incorporates the quartic Yukawa 
I d'^C'^v'^ of the oscillator coupling and the system-environment coupling of the form 
^gym^fX^- (We parametrized the previous /i by a dimensionless gy = fi/m^). 
Furthermore for simplicity we assume that the renormalized if mass vanishes. Since 
we take the lowest non-trivial order in perturbation theory, the rj dependence of A 
may be ignored in the discussion here. The inequality condition ( |5.29| ) is then, using 
(13:261) 

k gl 1 f Vk^ + 4m^ \ 

Physical processes that contribute to the dissipative friction here are X + ^ X 
that becomes dominant at high temperatures. If instead the dominant process of 
dissipation is — * X + X ? then the following constraint may be modified substan- 
tially. We hope to cover this case in subsequent investigation that involves non-linear 
dissipation. The growth rate A does depend on the mode variable, specifically the 
{k , g^o) value within the instability band. It was however shown in ref that for the 
large amplitude of g^o ^ f^s. ? the average A ~ 0.2 . The inequality thus limits the 
momentum value of produced particles. Since the quantity k/nin the left hand side 
at a fixed k decreases as the temperature T increases, the momentum range becomes 
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smaller as T increases. But ior k > T, k/n rapidly increases, as the momentum k 
increases. Thus this constraint is not very stringent for k >T. 

Thus 77 = 2Xm^ gives the critical friction below which the parametric amplifi- 
cation occurs, as illustrated in fig.6 — 8. Although the case of the overdamping 
1] > 2Xm^ is not our main concern, we shall write down results in this case, which 
have some interesting features worth to mention. Numerical results in this case are 
illustrated with the largest rj value (r] = 1.0) in fig.6 - fig. 8. At high temperatures 
the basic physical quantities behave asymptotically; 



S 



T 
2 



, 2T 

In — + 1 , 

V 



fo 



V) 



2ujqT 



(5.31) 
(5.32) 
(5.33) 



As expected, there is no effect of the parametric amplification. The S parameter 
that characterizes the spectral shape is 



V 



1 + 



^2 



(5.34) 



Time evolution of this spectral shape parameter is depicted in fig. 8 by = 1.0 line 
for the large friction. The oscillatory behavior around 6 = 0.9 is clearly observed. 
Thus even a large friction cannot destroy the quantum nature of the parametric 
amplifier: In a sense the system never completely equilibrates. 

We next turn to the low temperature limit of T ^ in which coth ^ — > 1 in the 
a;— integral. In the reheating problem this case is more relevant, since the universe 
supercools after inflation. In this discussion we need the high frequency cutoff at Q 
and flrst utilize with e = 1/ft , 

d s — T 

/o d{s — r) (s — r)^ + 



f 

Jo 



dijjijj cosuir — s) e 



(5.35) 



Singular terms may arise as e ^ 0+ . We must retain the leading logarithmic sin 



gularity such as ln(^^) 



ln(— ) , because this reflects accumulated effect to the 



high frequency cutoff f2. We shall explain estimate of the basic quantities only in the 
simpler case of the small friction, t] < 2Xm^, and state the flnal result for 1] > 2Xm^. 
Numerically both the small and the large friction cases have been computed. 
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Let us illustrate concretely some integration: 
A — — dr ds . . „ 

TT Jo Jo I'l + l^ 



u(t)u(s) r ^ , , (i S — T , 

^ ^ exp - 2i - T - s — . (5.36 



We partial integrate with respect to s and expand around s = r to estimate this 
type of integral. The leading term of A above is the surface term in the partial 
integration, which gives 

rt 111 t) , f — 7- 



1 [ 
TT Jo 



u{t) 



(5.37) 



The rest of the integral in A does not contain the logarithmic term in the e — > 
limit, and may be ignored. 
Similarly 

'■^ , v(t)v(s) r V ,^ N 1 fi- 

ds ^ ' exp[ - ^ (2i - r - s) ] — 

r 
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- fdr [ 

TT Jo Jo 
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TT Jo 
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dr ——^ e 

-r)t/2 



T2 + e2 
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d{s - r) (s - r)2 + e2 
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?7 u(t) 
2 ~ M(t) 



(5.38) 
(5.39) 

(5.40) 
(5.41) 

(5.42) 

(5.43) 

(5.44) 



The reason why the dependence on the initial state factor (3o disappears is dominance 
of the ln(;^) term in K. 

Physical quantities such as the average number of produced particles and the 
entropy are asymptotically 

1 T) iiit) 



uJo {N) 



E 
2" 



1 + 



u{t) 



,«(t)_^2 



'^ln(i^) e-'*/^|44 



(5.45) 
(5.46) 
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This result markedly differs from what we know without the environment effect: A 
pure quantum state has zero entropy, while interaction with environment gives a 
randomness in agreement with S ^ ^ ln{N) , for large (N). The energy distribution 
is given by 

1 p-W(2M) 

/M ^ -^==^^, (5.47) 

- 4ln(-) fl + A(^-44)0"- (5-48) 
In the large friction of ?7 > 2Am^ at low temperatures a similar analysis gives 

{uj) ^ ^My), (5.49) 
S = S(x), x^-J81n(— ), (5.50) 



with S'(x) defined by eq.( [4.67|) . 

The formulas written here are valid for t] In(^) ^ luq. When this does not hold, 
the factor, r] In(^) , should be replaced by some intrinsic frequency scale of the 

system. Our results at low temperatures indicate that the difference from the high 
temperature case is described by a replacement, — > In(^) , ignoring subleading 
terms. Thus the temperature plays a minor role in the parametric amplification. The 
crucial condition is the inequality, 2X{ri) > r] . 

Results of numerical computation at T = are shown in fig. 11 — fig. 13. Global 
features of time evolution at low temperatures do not differ much from the high 
temperature case, as anticipated by analytic results above. Thus if time is displaced, 
most physical quantities at different temperatures look the same. This is a surpris- 
ing, and an important result of our work, indicating that the parametric effect is 
insensitive to the environment temperature except through the rate 2X{ri) — rj . 
In fine details this universal feature is somewhat spoiled since the amount of the 
time displacement is not the same for the average energy (u) and the entropy 5*. 

As advocated in another model |p , our results indicate that the energy spectrum 
at the instant of particle creation deviates from the thermal one: it is neither of 
the form of Bose condensation. However, the spectral shape does not necessarily 
persist till the last moment, since interaction among created particles neglected in 
the present work can establish thermalization, as discussed in ref 0. 
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Our computation in the present work has been hmited to a single mode by spec- 
ifying k and (the amphtude of ^ oscillation). To assess more fully effect of the 
parametric amplification, one must further mode-sum over k in all possible instability 
bands. A consequence of this mode-sum might be to erase the oscillatory behavior in 
time evolution of the average energy and the entropy. Thus although the quantum 
behavior caused by the parametric amplifier persists in medium, the mode-sum may 
or may not erase the bizarre oscillatory behavior observed in a single mode calcula- 
tion here. We have also performed a sample computation when the parameters of 
the coherent oscillation are outside the instability band and inside a stability region. 
The exponential growth is not observed in this expected, but the system does 

not amalgamate to the environment, either. The quantum behavior such as a per- 
sistent oscillatory behavior is observed, thus it seems that the periodic perturbation 
keeps the quantum system from losing its quantum nature even in thermal bath. 

Although it is not our main topic to deal with the reheating problem after in- 
flation, it might be of some interest to comment on this. Our approach is based 
on the premise that separation between the system and the environment is evident. 
In the preheating stage after inflation, especially prior to copious production of if 
particles, the separation is obvious. Indeed, the relevant environment is the vacuum 
at zero temperature. Whether the Bose condensation may or may not occur can be 
decided at this stage. On the other hand, after explosive particle production the 
environment changes. Clearly the new environment consists of created particles due 
to the infiaton decay prior to thermalization. We cannot discuss this stage in our 
framework. This is regrettable since this stage is the most important to the later 
cosmological evolution. We also note that we have neither discussed how the infia- 
ton decay is completed: here the back reaction caused by particle production is very 
important. Some recent works on particle interaction after their creation are found 

n, n, n. 
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Although our achievement on the reheating problem is very modest, our result 
on the environment effect to the parametric amplification may be of considerable 
use in other cosmological problems. These are the Polonyi or the moduli problem. 
These weakly interacting fields are thought to decay after long periods of matter 
dominance, thus potentially giving a serious entropy crisis at their decay pSf. Their 
oscillation lasts long enough in cosmic medium, thus it is important to estimate the 
environment effect. 
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6 Summary 



Since a great many have been described in the present paper, it would be useful to 
summarize here what we think we achieved in this work. Our work addresses how to 
quantify the effect of the environment on the behavior of a quantum system and can 
be classified into three parts: The first concerns general formalism to characterize 
and approximate the system-environment interaction; the second concerns the effect 
of the environment on the general dynamical system of time dependent harmonic 
oscillators; the third concerns the effect of the environment on dynamics of harmonic 
oscillator with periodic frequency. 

Our results on these three issues are as follows. 
On general formalism of the system-environment interaction: 

• We extended the influence functional method by introducing the spectral weight 
r{u}), called the response weight. For the bilinear coupling between the sys- 
tem and the environment variables the correlation function that appears in the 
influence functional is identified to the Green's function in the environment 
medium, for instance the real-time thermal Green's function for the environ- 
ment in thermal equilibrium. The response weight completely characterizes 
the system-environment interaction. 

• Wc developed two methods of modeling the system-environment interaction. 
The first one is a mcromorphic approximation to the response weight, some- 
what akin to the resonance approximation in the low energy scattering problem. 
The second one is a perturbation theory in terms of the interaction strength. 
The Ohmic behavior at low frequencies, r{u}) oc a;, is a common feature of 
models discussed in the present work. 

• We introduced a class of non-Gaussian field theory models of the interaction 

Lagrangian density of the form £int = — | l^ffX^ between the system field (p and 
the environment field x- Finite temperature computation in this field theory 
model has been worked out at one loop level. How to model this case as an 
extended Ohmic model was briefiy discussed. 

On the time dependent harmonic oscillator in thermal bath: 
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The late time behavior of the system under the environment effect is described 
by an approximation to the correlation function; the dissipative part of the 
correlation given by a local form, ^/(r) ~ ^?5'(t) , where the friction constant 
rj describes the inverse time scale of relaxation of the system variable towards 
thermalization. Correspondingly the response weight r{u!) = ^uj f{^) , with 
f{x) decreasing rapidly as x — * oo. Q is the cutoff of the high frequency part. 
We identified the response weight and the friction constant in the two models 
of the system-environment interaction. 

We solved in a closed form the problem of how time dependent harmonic 
oscillators behave asymptotically at late times: The reduced density matrix of 
the system under the effect of thermal environment is determined in terms of 
two quantities, a single classical solution p(t) obeying P —rjp + u!^{t)p — 
with p(0) = , p(0) = 1 , and the response weight. 

Physical quantities such as the energy distribution, the entropy of the sys- 
tem, and the total population number in oscillator levels have been explic- 
itly computed. The resulting energy spectrum is given by a universal func- 
tion of the scaled energy, x — energy / (average energy), having the form, 
exp[—x/S^] Iq{x \/Y^^^ / 5"^ ) / 5 , where the time dependent parameter 1 — 5 
(with < 5 < 1) describes a degree of departure from thermal equilibrium of 
the system. 

the environment effect against the parametric oscillator: 

For the first time it was verified that the parametric amplification can occur in 
medium, using the infiuence functional method. The average energy and the 
total number of population exponentially increases with time, if the growth 
rate of the parametric amplification in medium is larger than the friction 77, 
interpreted to be a measure of the relaxation rate towards thermalization of 
the system variable. The critical friction for the exponential growth has also 
been given, above which the parametric effect does not occur. 

The resulting energy distribution deviates from the thermal Boltzmann dis- 
tribution, having a long tail of the high energy component characterized by 
an exponentially increasing average energy and an exponentially decreasing 5 
parameter. 
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• Late time behavior of the average energy, the entropy, and the energy spectrum 
is insensitive to the environment temperature if time is displaced. This is true 
even including the T = vacuum. 

After submission of this paper for publication, we became aware of some related 
works on the reheating and the preheating problems after inflation pGf , ||2^ , . 



47 



Acknowledgment 

This work has been supported in part by the Grand-in-Aid for Science Research 
from the Ministry of Education, Science and Culture of Japan, No. 08640341 and 
08740185. One of us (M.H.) acknowledges the hospitality of the SLAC theory group 
while part of this work has been done during his stay there. 



48 



Appendix A How to numerically compute A,B ,C 



The basic quantities we must know before computing physical quantities are A,B ,C 
introduced in Section IV. These contain some factors hke u{t) in denominators, and 
although zeros of these functions cancel in the end, it would be much better to avoid 
these zeros in intermediate steps of calculation. We achieve this by introducing, and 
expressing A,B ,C in terms of, supplementary functions that obey simple differential 
equations with definite boundary conditions. We briefly describe here the method 
to do this. 

Let us first recall that the fundamental function u{t) obeys the following rela- 
tions; 

u{t) = f{T)u{T), (A.l) 

u(0) = 0, u(0)^l. (A.2) 

For notational simplicity we replaced the time dependent frequency uj'^{t) by — /(r). 
The secondary function v{t -,1) is then 

v{t;t) = 0, v{0;t)^-l. (A.4) 

We distinguish the two arguments, r and t, to indicate that r is a time variable and 
t is the end time at which we set the boundary condition. The dot always means 
the r derivative. The normalization here was chosen such that the Wronskian 

u{T)v{T]t) -u{T)v{T;t) ^1. (A.5) 

We further define two more functions; 

k(T-t) ^ «(f) ( - i.(T) + I ) 

Note that t > r. These functions are shown to satisfy 

d^h{T;t) = f{t)h{r-t), (A.7) 
M0;0) = -1, aiMO;0) = |. (A.8) 
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Finally we define 



cit) 

s{t) 

which are shown to obey 



= - u{t) f dTv{r; t) cos{ut) e"^^^ , 
Jo 



— u{t) / dr v(t ■,t) sm(ujT) e 
Jo 

lim hlr ; t) , 



7jr/2 



9(0) 

dMt) 

c(0) 

s{0) 



1 , dtg{0) 



2 ' 



f{t) c{t) + cos (cut) 6""/^ , 
0, dtc{0)^0, 
f{t) s{t) + sm{ut) , 
0, dts{0)^0. 



(A.9) 

(A.IO) 
(A.ll) 



(A.12) 

(A.13) 
(A.14) 
(A.15) 
(A.16) 
(A.17) 



We introduce the following multiplication and the integral transform to simplify 
subsequent expressions; 



A = 



Vo = 
Basic functions are then 



yCOth(^) X , 

°° 3ui 

dujr{uj) coth( — ) x F{uj) 
2 



s[Vo{u\t) + S^)+V^{c\t) + s\t) 



(A.18) 
(A.19) 

(A.20) 



^+VoV^ ((eit) + ^^)ic\t) + s'{t)) 



+ ( ^it) + ^ ){c'{t) + s^t) ) - 2( u{t)u{t) + 



){c{t)c{t) + s{t)4t)) 



+ {V^{ c'{t) + s'{t) ( c\t) + s\t) )-iVU cmt) + smt)) Y} ] , (A.21) 

C = AAe-^ (vo{u{t)u{t) + ^^^^)+V^{c{t)c{t) + m<t))]-'l. (A.22) 
V '^o / 4 

We thus need to solve four linear differential equations for u{t) , g(t) , c{t) , s{t) , 
that obey different boundary conditions, to get A,B ,C. Some of these equations 
contain inhomogeneous terms that depend on cu. 
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Figure captions 



Fig.l 

Response weight in tlie field tlieory model at finite temperature. One loop result 
of r{uj) = 2QIl{uj) with S5n(ct;) the imaginary part of the self-energy diagram is 
shown for a parameter set; T = 50,k = 5,fi = l in the energy unit of m (the mass 
of decaying particle). Also shown by the dotted line is the approximation to this 
weight given by eq.(3. p.34D in the text. 



Fig.2 

Energy distribution given by the universal form, /o(x ; S) of eq.(4. [4.56| ), plotted 
against a scaled energy of x = uj/ {uj) for three choices of 5 = 1.0 (the dotted line), 
0.5 (the solid line), 0.2 (the broken line). 

Fig.3 

Time evolution of populated number (= average energy /wq — |) of harmonic 
oscillator levels for various choices of the parameters, r] (the friction) and Tq (the 
system temperature), all at the environment temperature T = 100; 
iv , To) = (0.1 , 0) , (0.2 , 0) , (1.0 , 0) , (0.2 , 100) , (0.2 , 150) 

in the energy unit of Uq, given by the designated lines with a, b, c, d, e, respectively. 
Fig.4 

Time evolution of entropy for harmonic oscillator for the same set of parameters 
as in Fig.3. 

Fig.5 

Time evolution of the fundamental solutions for the modified Mathieu equation, 
normalized at r = t to give (the solid line), and at r = to give (the dotted 
line). Relevant parameters are h = 1 , 26 = 0.75 , rj = 0.2 , t = , and 

the time scale is measured by 

Fig.6 

Time evolution (in the same time scale of Fig.5) of occupation number (= av- 
erage energy /uq — ^) of oscillator levels of reference frequency Uq in the case of 
the parametric oscillator with h = 1, 26 = 0.75. Three different situations are 
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compared: 77 = or in vacuum by the broken line, the small friction case with 
{r] ,T) = (0.2 , 100) by the solid line, and the large friction case with (1.0 , 100) by 
the dotted line. The system temperature is all set at Tq = 0, the ground state. 

Fig.7 

Time evolution (in the same time scale of Fig. 5) of entropy in the case of the 
parametric oscillator with h — 1 , 29 — 0.75 . Three different situations of the 
friction rj and the environment temperature T are compared: {1] ,T) = (0.05,100) 
by the broken hne, (0.2 , 100) by the solid hne, and (1.0 , 100) by the dotted line. 
The system temperature is all set at Tq = 0, the ground state. 

Fig.8 

Time evolution of the spectral shape parameter 6. Two cases of the small and the 
large friction at a high temperature (T = 100) for the parametric oscillator {h — 1.0 
and 2 9 — 0.75) are compared to the simple harmonic oscillator (with 77 = 0.2) 
depicted by the dotted line; the small friction case (77 = 0.2) given by the broken 
line, and the large friction case {rj = 1.0) given by the solid line. 

Fig.9 

Time evolution of the basic quantities A, B in figure 9A and C in figure 9B, 
for /i = 1 , 2 e = 0.75 , 77 = 0.2 , T = 100 . The time scale is measured by the 
angular oscillation frequency of the parametric oscillator, 

Fig.lO 

Time evolution (in the same time scale of Fig. 5) of integrated energy fiow from 
the parametric oscillator into the total system including the environment. Three 
cases are compared: 77 = (isolated system) by the sofid line, (77 , T) = (0.2 , 100) 
by the broken line, and (1.0 , 100) by the dotted fine, all at the system temperature 
To = 0. Inside the inlet is shown the case of the large friction of (1.0 , 100) in a linear 
enlarged scale. 

Fig. 11 

Time evolution (in the same time scale of Fig. 5) of populated number of oscillator 
levels of a reference frequency cuq in the case of the parametric oscillator with h = 
1 , 29 = 0.75 . Three different situations in vacuum (Tq = 0) are compared to the 
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high temperature case by designated letters: {rj , T , Tq) — (a) (0.2 , , 0) ; (b) 
(0.2, 100, 0) ; (c) (0.2, 0, 100) ; (d) (1.0, 0, 0) ; (e) (0, 0, 0). 

Fig.l2 

Time evolution (in the same time scale of Fig. 5) of entropy for the set of param- 
eters; (?7 , T , To) = (a) (0.2 , , 0) ; (b) (0.2 , 100 , 0) ; (c) (0.2 , , 100) ; (d) 
(1.0, 0, 0). 

Fig. 13 

Time evolution (in the same time scale of Fig. 5) of the spectral shape parameter 
S with the friction rj — 0.2, at high (T = 100 by the broken hne) and at low (T = by 
the sohd line) temperatures, taken for the parametric oscillator oi h— 1, 29 — 0.75 . 
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